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<N ■ Abstract 

Sphaleron energies within the standard model with a real Higgs singlet added on 
are calculated. The coupled non-linear equations of motion are numerically solved 



and the sphaleron energy evaluated for a set of parameters in the Higgs potential. I 
find a small difference in the sphaleron energy compared to the standard model. A 
slightly stronger constraint on the strength of the first order phase transition thus 
results for this model. 

1 Introduction 



X 

If the baryon asymmetry of the universe was generated at the electroweak phase transition, 



as many models are now speculating [JTJ] , it is important that the rate of baryon violation be 
determined more accurately. This rate in the broken symmetry phase is directly related to 
the energy of the sphaleron configuration in the model, because the sphaleron energy sets 
the height of the energy barrier between the topologically inequivalent vacua, and baryon 
violation occurs whenever such a vacuum transition takes place. To preserve any baryon 
asymmetry created at the phase transition the baryon violation rate must be suppressed 
in the broken phase and this imposes a constraint on the sphaleron energy at the phase 
transition. Hence it is important to calculate the sphaleron energy within each model 
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being considered. It is in any case interesting to know the explicit form of the sphaleron 
configuration in extensions of the Standard Model. 

Manton et.al calculated the sphaleron energy within the Standard Model (SM) using 
a radially symmetric ansatz for the gauge and Higgs fields 0. They obtained 

E sph = —B(\/g 2 ) (1) 

where v is the Higgs vacuum expectation value (VEV), g is the weak gauge coupling and A 
is the quartic self coupling of the Higgs boson. Their calculations showed 1.5 < B < 2.7 for 
X/g 2 between and oo. Since the baryon violation rate in the broken phase is proportional 
to e~ Esph ' T , suppression of the baryon violation rate in the broken phase requires B 



E{T C ) > 



T, 



45, (2) 



c 



independent of the model, where T c is the critical temperature for the phase transition. 
Hence in the SM, one obtains 

£ ~ 1-4 (3) 

J- c 

as the condition on the strength of the first order phase transition required. 

Because this constraint seems so restrictive in the SM [apart from other possible in- 
adequacies like the strength of Kobayashi-Maskawa (KM) CP violation] it seems likely 
that an extension of the SM is required to generate sufficient baryon asymmetry at the 
electroweak phase transition. Hence it is important that the sphaleron energy be calcu- 
lated in these models. In this paper I consider the sphaleron solutions in the minimal SM 
augumented by a real singlet Higgs field. This model was shown to be capable of satisfy- 
ing at least one of the several conditions needed for electroweak baryogenesis, namely the 
condition that any baryon asymmetry created at a first order phase transition does not 
get washed away by subsequent sphaleron processes ||. This model is interesting due to 
the existence of tree-level trilinear Higgs couplings which seem to enhance the possibility 
of satisfying this condition (I also point out other possible interest this model has in this 
regard in the discussions later). Kastening and Zhang || have investigated sphaleron 
solutions when a complex Higgs singlet is added to the SM (their paper also applies to 
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the real Higgs singlet case but with the trilinear coupling terms missing) and found non- 
zero differences in E sp h compared to the SM case. I too want to compare the sphaleron 
energy of this model with the SM. If the calculations here show a marked increase in the 
sphaleron energy compared to the SM value then one can expect a greater suppression 
of the sphaleron rate in the broken phase, relaxing the constraints on the Higgs masses 
somewhat. If on the other hand a considerable decrease in the sphaleron energy is seen 
in this model, then the first order phase transition will need to be even stronger and 
consequently less parameter space will be available to prevent a washout of the baryon 
asymmetry of the universe. The present work is thus a refinement of Ref.0], since the 
minimal SM sphaleron configuration was used there as an approximation. 



2 The equations of motion 

A real singlet Higgs field S is added to the minimal SM with its doublet Higgs field 0. 
The Lagrangian of the gauge and the Higgs sector of the model is then: 

C = —F^Ff + (D^)\D»<P) + l -d»Sd»S - V(<p, S), (4) 

where 

K v = d"W v a - d v W£ + ge ahc W£W v c (5) 

v(<p, s) = a,(0W - + y s * - Y s2 ~ f s " + 2X ^<f>)$ 2 ~ \ WW. (6) 

= (d, - UgT a Wpcf>. (7) 

(For simplicity, the fermion fields will be set to zero as well as the Uy(l) gauge field by 
setting the weak mixing angle to zero.) 

The Higgs potential V(<f>, S) breaks the SUl(2) <g> Uy(l) symmetry to Uq(1) when the 
Higgs fields take on the VEV's at: 

/ 



, and (,(,)) = 
u I V ^ 
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where u = 246 GeV and v is kept non-zero for full generality. The minimum of the 
potential (pi) is given by (u, v) where u, v are solutions to 



(-/ij + X^u 2 + Xv 2 - av/2)u = 

— [igV + Xsv 3 + Xvu 2 — au 2 /A — av 2 = 0. 

I eliminate the parameters //| and /x| in favour of one of the (u,v) solutions to Eq.( 



(9) 



^ = X^u 2 + Xv 2 - av/2 

fi 2 s = Xsv 2 + Xu 2 — au 2 /Av — av. 



Now for static classical fields, one has 



(10) 



d n Fr = <9 O = d S = and W£ = 0. 



'11} 



Then the equations of motion derived from L are: 



2 

[2A (0t 
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+ 2AS'0 - -50 



[8A0 f - 2/4] S + 4A S S 3 - (70V - 2aS 2 . 



(12) 

(13) 
(14) 



And the corresponding energy functional is: 

£ = | cPx[±F%F% + (A0) f (A0) + ^X^) + 5)]- 

Now I use the spherically symmetric ansatz of the form: 

2i 



where 



W?T a dx l 



u 



-figV^dU 00 ^ 00 )- 1 
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( 



— h(gVr)U c 




1 



u 



X - 1 

r 



s = 

z x + iy 
—x + iy z 



V2 



p(gVr), 



and 



v = ^FT 



(15) 



(16) 
(17) 



Then the equations of motion become: 
= 2/(1 -f)(l- 2f) 





? de 






'?% 


— 


d 




di 



ev 2 

Av 2 



h\l-f) 



2/i(l - f? + (v 2 /* 3 - Ah + XvYh - ^=ph 



(19) 



4£ 2 



A, 



av 



V 



au 



-h 2 



where £ = gVr. And the energy functional becomes 



E = *f f °° di { 4 (|) 2 + _ /)]2 + ^ (f ) 2 + - /)]' 



+ 



_i_ A s t) 4 4 _ ra 2 n i|2 
2 P V2 P 



x^h 2 - /4/a^) 2 + (2A M 2 /i 2 - /4>y 

Co} , 



(20) 



fty^c 3 „3 
3 ^ 



where e Q is the normalisation needed to set the minimum value of V(<f>, S) at zero. 

Hence within the validity of the ansatz being used, the problem of finding the sphaleron 
energy in this model becomes equivalent to either minimizing the above energy functional 
( |2T)| ) or solving the coupled non-linear differential equations (|19D . The latter method will 
be used in this paper, by numerically solving the equations and then calculating the value 
of E which corresponds to the solutions found. 

Boundary conditions for the functions are required and to obtain these at £ = 0, the 
equations are Taylor expanded about £ = and regularity is imposed there to obtain 



£^0: f^a s £ 2 , h->b s £, p -> c s + c 2 — 



(21) 



where 



c 2 



— ( -j=u 2 + A\u 2 v - V2av 2 {l - c s ) + 2X s v 3 {l - c 2 ) 
v2 



(22) 



3g 2 vV 2 

while at £ = oo the Higgs fields are made to take on their zero temperature VEV's by 
f,h,p — > 1. The equations are actually singular at both boundaries, so for large £ one 
looks for asymptotic forms of the functions. I find 



oo : 



/ - 1 - a^, h - 1 - &,e-"7£, P - 1 - Qe- c V£ 



(23) 
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where 



' u 

w 2 ' 
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V2 2V2v 



(24) 



V g 2 V 2 ' 

There are three second order ODEs, so one needs to have six boundary conditions. The 
variables (a s , b s ,c s ,ai,bi,ci) implicitly contain these boundary conditions so the solutions 
to the equations of motion and these variables need to be determined simultaneously. 



3 Solutions and results 

Because there is no analytical way of finding the solutions to the coupled non-linear 
ODEs of (|i~9D, numerical methods must be used. In order to do this, some values of the 
parameters in the potential need to be chosen, and I choose the same set here as was 
used in Ref. for easy comparison. These sets of parameters then indicate the values 
for which the first order phase transition is strong enough while the mass of the smallest 
Higgs is above the experimental lower limit, where the minimal SM sphaleron energy is 
used as an approximation. Here I improve on this by using a corrected sphaleron energy. 

I used the shooting method to solve the equations of motion. The results obtained for 
each set of parameters are tabulated in Table 1, and a typical solution is shown in Fig.l. 



Note that the values of E are given here in units of AitV/g where V = y/u 2 + v 2 , unlike 
the SM values which are measured in units of 4nu/g. The results show that the sphaleron 
energies in this model are somewhat lower than the SM CclSG clS calculated by Manton 
et.al 0. Similarly when compared to the SM + a complex Higgs singlet case considered 
by Kastening and Zhang ||, I find lower values of the sphaleron energy. 

What does this mean for electroweak baryogenesis? Satisfying Eq.(0) in this model 
requires 

t c ~^b-' (25) 

where B is now a function of the parameters in this model. Table II compares the values of 
Abg/AirB obtained here with the V/T c values obtained in Ref 0] at the phase transition, 
for each corresponding set of parameters. They show that the condition (p5|) is easily 
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still satisfied, because the first order phase transition is still too strong compared to the 
decrease in the energy of the sphaleron found. 

Of course one must keep in mind all the uncertainties involved in such an analysis. 
The rate of sphaleron transitions in the broken phase involves prefactors to the Boltzmann 
factor || which is difficult to calculate near the transition temperature due to infrared 
divergences. I have also extrapolated my zero temperature sphaleron energy up to a 
high temperature in making the comparison ( J25|) M and have not taken any quantum 
corrections into account 0. Also, I have not considered any source of CP violation in 
this model at all. There is still some debate on the question of whether or not the SM 
can provide sufficiently strong CP violation in order to produce the observed baryon 
asymmetry ||. If it turns out that it can't, this model being considered here will need a 
new source of CP violation || . If however it turns out that the SM can provide sufficiently 
strong CP violation, then adding the real Higgs singlet will be an extremely economical 
way of achieving electroweak baryogenesis as it can provide a first order phase transition 
for a larger region of parameter space than can the minimal SM (Once again there is much 
uncertainty in the calculations of the effective potential at finite temperature so hasty 
conclusions cannot be made.) It also does not suffer from the possible problem faced by the 
two Higgs doublet model which may not be able to produce sufficient baryon asymmetry 
in the first place even though it can prevent a washout of any baryon asymmetry that 
has been created [jlTJ. This is because there are two VEV's and the phase transition can 
be in two stages. If after acquiring the first VEV the sphaleron rate is too suppressed 
because the gauge symmetry is broken, not enough baryon asymmetry will be created 
before the second phase transition which will preserve this asymmetry. Even though the 
real Higgs singlet model considered in this paper also has a two stage phase transition, 
because the gauge symmetry remains unbroken after the singlet acquires its VEV, the 
baryon violating rate will not be suppressed in between the two phases. 
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Table 1: Representative parameter values and the corresponding sphaleron energies (in 
units of AnV/ g). The parameters A^s and A are dimensionless while the other parameters 
have dimension GeV. 





As 


A 


V 


a 


a 




0.071 


0.082 


-0.0234 


-125 


-12.81 


12.20 


1.40 


0.063 


0.073 


0.0250 


-154 


-13.26 


3.58 


2.20 


0.087 


0.090 


0.0114 


-245 


-9.37 


1.24 


0.82 


0.096 


0.082 


0.0499 


-146 


-17.27 


5.10 


1.80 


0.050 


0.077 


0.0168 


-117 


-12.43 


8.34 


1.91 


0.057 


0.0920 


0.0026 


-171 


-16.39 


4.64 


0.55 


0.057 


0.0920 


0.0026 


-107 


-16.54 


16.76 


0.91 



4 Conclusion 

I calculated the sphaleron energies in the standard model with a real Higgs singlet added 
on, for some representative sets of parameters. These show a decrease in the actual 
sphaleron energies compared to the values found in the Standard Model. However the 
decrease is not sufficient to constrain the strength of the possible first order phase tran- 
sition in this model. If KM CP violation in the Standard Model is sufficiently strong for 
electroweak baryogenesis, then this model will be an extremely economical way of creating 
the baryon asymmetry of the universe at the electroweak phase transition. 
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Table 2: Comparison of V/T c obtained from constraint (2) and E^ found with V/T c 
obtained at the first order electroweak phase transition. 



45g/4nB m^ s 



V c /T c 



4.10 



2.75 



2.48 



1.25 



1.18 



1.03 



1.61 



103, 108 



61,102 



75, 109 



67, 79 
73,95 
68,92 



73,88 



170/107 
390/93 
510/117 
430/97 
360/70 
360/85 
360/63 
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Figure captions 

Fig.l: A typical solution to the coupled non-linear ODE's of Eq.(19). (The fourth line of 
parameters in Table 1 was used here.) 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409360v2 




gvr 



